Probing neutrino mixing angles with ultrahigh energy neutrino telescopes 
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We point out that detecting v^s hom distant astrophysical sources with the up-coming and future 
neutrino telescopes using the Glashow resonance channel Vc.e~ W~ — > anything, which occurs over 
a small energy window around the energy of ~ 6.3 PeV, offers a new way of measuring or setting 
limits on neutrino mixing angles, in particular the angle 9i2, thereby providing an independent 
experimental probe of neutrino mixing angles. We also discuss how this exercise may throw light 
on the nature of the neutrino production mechanism in individual astrophysical sources. 
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Experimental determinations of the parameters gov- 
erning neutrino flavor oscillation phenomena, such as the 
mass-square differences and mixing angles amongst the 
three light neutrino species, have so far been done with 
solar- atmospheric-, reactor-, and accelerator neutrino 
experiments [l| that involve neutrinos of few MeV - GeV 
energies. At a purely phcnomenological level, one may 
ask if neutrinos of vastly different energies oscillate in 
the same way. In this context it is important to think of 
other possible independent ways of measuring the neu- 
trino parameters that involve neutrinos of different (in 
particular higher) energies than those currently employed 
in neutrino oscillation experiments. 

In this Letter we point out that up-coming and future 
neutrino telescopes capable of detecting neutrinos of very 
high (at least up to several PeV: 1 PeV = 10^ GeV) en- 
ergies, would potentially also be able to measure one or 
more of the neutrino mixing angles at such high ener- 
gies 0. Specifically, we suggest that detecting j>e's from 
distant astrophysical sources through the Glashow reso- 
nance (GR) 3\ channel, i'ee~ — > W~ — » anything, which 
occurs over a small energy window around the De energy 
= m^/2'me = 6.3 PeV, offers a new way of measur- 
ing the angle 9i2 (see below) 0. 

Neutrino flavor oscillation phenomena are governed by 
six independent parameters 5J: two mass-squared differ- 
ences, Amfj = — m\\ and Am^^ = | m| — m^l, 
three mixing angles, 612 , 6*23 , and 6*13, and a possible CP- 
violating phase S. Here (mi , m2 , rus) are the masses cor- 
responding to the three light neutrino mass eigenstates 
(i/i ,V2 ,1^3). The mixing angles parametrize the relation 
between the mass eigenstates and the flavor eigenstates 
(ue ,v^,Vt)', see equation © below. The current exper- 
imental situation, as summarized in Ref. 0, is the fol- 
lowing: the solar neutrino data are consistent with flavor 
oscillations mainly between v^. and driven by a mass- 
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mixing angle ^soiar = ^12 ~ 32°. 5, while the atmospheric 
neutrino data are explained by oscillations mainly be- 



tween and Vr with Am^+ 



Am^3 



2.6 X 10--^ eV 



and 6'atm = ^23 ~ 45°, the maximal value. The other 
mixing angle ^13 is constrained by reactor neutrino ex- 
periments to be very small, 6113 <, 9°. The CP-violating 
phase 5 remains undetermined by present experiments. 
The uncertainties (currently roughly 10-20%) in the val- 
ues of the above parameters may be expected to be sig- 
niflcantly reduced by future experiments. 

Very high energy neutrinos of energies well above 
several hundreds of TeV, and in some cases extend- 
ing well into EeV (10^ GeV) energies, are predicted 
to be produced by astrophysical sources such as Ac- 
tive Galactic Nuclei (AGN) and Gamma Ray Bursts 
(GRB)'s 0. More "exotic" sources such as cosmic 
topological defects are predicted to produce neutrinos 
of energy even up to several hundred EeV 0. De- 
tection of such high energy neutrinos from distant as- 
trophysical sources is the primary goal of a host of 
proposed or already-under-construction large under-ice 
and underwater optical Cerenkov detectors such as ICE- 
CUBE 8], ANTARES 0, NESTOR 0, NEMO Jli, as 
well as radio Cerenkov detectors such as RICE 12] and 
ANITA jl3|. These detectors wiU be able to determine 
the energy as well the arrival directions of the neutrinos 
enabling point source neutrino astronomy. 

High energy neutrinos are produced in astrophysical 
sources mainly from the decays of charged pions (and 
kaons) which are produced through interaction of high 
energy protons with either ambient photons ( "7)7" ) or 
protons ("pp") within the source 14]. In the P7 process, 
charged pion (and hence neutrino) production occurs 
dominantly through production of the A resonance Il5j: 
P7 — > A+ TT+n, with 7r+ — s- /i^f^, and /i+ — > e+i^e//^. 
Thus each p7 interaction eventually produces one each of 
i/g, v^, and Vf^. Denoting by Ty = {v^ ,Ve ,Vr .Vt} 

the fractional amounts of j/'s and i?'s of various flavors, 
we thus have for the P7 process the flavor fractions at 



the source, J^p^ = {4,0 



0,0}. In the pp process, 



on the other hand, each inelastic pp collision produces a 
nearly equal mix of tt+'s and 7r~'s whose decays together 
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produce the flavor fractions, T^p = {^,^,^,^,0,0}. 
Note that the (v+v) fractions, J-(^^^^j = {{ve + Ve) , + 
Dfj^) , {i^r + i^r)} — {| 1 1 1 0} for both pj and pp processes. 

In the tenuous (radiation dominated) environments 
within most potential astrophysical sources such as 
AGNs and GRBs, the dominant production mechanism 
for high energy neutrinos is expected to be the pj pro- 
cess However, for generahty let us assume that a 
fraction x of the total number of all v's and P's produced 
at the source are of pp origin and the rest of p^ origin. 
Thus, a; = for a pure p'j origin of the neutrinos. The 
flavor fractions at the source now work out to be 
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(1) 



Note that (y + v) fractions of different flavors at source 
are again T^i^^^ = {■! , | , 0}, independent of the fraction 

X. 

The flavor fractions of the i^'s and j>'s arriving at 
Earth are different from the above fractions at the 
source due to neutrino flavor oscillations during the 



propagation from the source to the Earth. In the 
standard, 3-generation scenario of neutrino oscillations, 
the flavor eigenstates \va) (a = e,/i,r) are related 
to the mass eigenstates \vi) (i = 1,2,3) through 
\va) = J2iUaiWi)- For values of Am^ as indi- 
cated by solar and atmospheric neutrino experiments, 
the vacuum oscillation length Lose ~ AnEi,/\/^m?\ ~ 
2.5 X 10-12 {E^/lTeY) (lO^^ eVVlAm^l) Mpc is always 
much smaller than the distances (several hundreds to 
thousands of Mpc) to possible astrophysical sources of 
neutrinos such as AGNs and GRBs, even for the highest 
energies of interest. The neutrinos thus oscillate many 
times before reaching Earth, and the oscillation (or sur- 
vival) probabilities averaged over many oscillations take 
the simple form 



P{^c ^^13)= I I' I I' ■ (2) 

i=l,2,3 

A "standard" form |l'| of the neutrino mixing matrix U 
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. Uel Ue2 Ue3 \ / C12C13 S12C13 8136 

U= \ Uf,i Uf,2 Uf,3 \ = \ -S12C23 - ci2S23Si3e'* C12C23 - si2S23Si3e**' S23C13 X diag [e'"^'/'^ , e'^'^''/^ , I 
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where Sy = sin6'y , — cosfly , and 5 ,ifi, and (^2 are 
CP-violating phases. The phases ^pi and Lp2 are present 
only for Majorana neutrinos; they do not enter into the 
expressions for oscillation probabilities [3, although the 
"Dirac" phase 5, if it is non-zero, does. 

In our numerical calculations below we use the above 
full form of the mixing matrix (|3J) in calculating the 
expressions for various oscillation probabilities |(5J with 
5 set to zero for simplicity. However, simple analyti- 
cal derivation of our main results are possible if we set 
9i3 ~ 0° and 6'23 = 45° (conditions which we hereafter 
refer to as "optimal"), in which case the mixing matrix 
U takes the simple form 
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Using Equations (^3), ©j and Q we can write the 
flavor fractions of i^'s and i^'s arriving at Earth in the 
optimal case as 

TrEarth ^ / 1 _ ^ ^ + g ^ -g 1 + g l_g 1 

-/-..optimal \l^3 6 ' 6 ' 6 ' 6 J 

(5) 



where p — i (x -I- 2^ — 2x^) , q = — xS^) , with ^ = 
s?2C?2 = isin^ 26*12. 

Note that (j^ -I- D) fractions of different flavors at Earth 
are simply {i , i , i}, independent of the fraction x in 
the optimal case. Also, the {v + v) fractions are indepen- 
dent of the neutrino mixing angles in the optimal case. 
Dependence on x and mixing angles (only 9i2 in the op- 
timal case) is contained in the separate number fractions 
of i^'s and v^s of individual flavors. However, separate 
identification of the i/'s and v^s of individual flavors is 
not generally possible through the usual charged-current 
(CC) interactions of neutrinos in the currently operating 
or planned water or ice-based detectors, since the sign of 
the charge of the produced charged lepton cannot be de- 
termined. The only exception is the which can be iden- 
tified in water or ice based detectors through the Glashow 
Resonance (GR) channel already mentioned An ex- 
perimentally measurable ratio such as 9^/ {Vfi + Vfj) giving 
the i>f. fraction to (z^^ +i^p) fraction measured over an en- 
ergy interval centered at the GR energy can then yield a 
measurement of the angle 9i2. Note that taking the ra- 
tio as above cancels out the unknown total neutrino flux 
coming from the source. Alternatively, if the spectrum 
of the neutrinos is the same across a broad range of ener- 
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Fraction of all neutrinos from pp interactions Fraction of all neutrinos from pp interactions 



FIG. 1: The ratio of number of Glashow resonance events 
to {i/fi + v^) charged current events in a water or ice based 
detector as a function of the fraction x oi {y -\- i/)'s that origi- 
nate from pp interactions at the source. The sohd hues, from 
bottom to top, are for 612 = 0° , 15° , 32.5° , and 45° — all 
for the "optimal" case 623 ~ 45° and 613 — 0°. The dashed 
lines are in the same order of 612 values as the solid lines — 
all for 013 = 0°, but 6*23 = 40°. 



gies, then one can also get the absolute normalization of 
the total neutrino flux from the total number of (j/^ + D^) 
events over a suitably chosen range of energies. 

The cross section for the GR interaction of De can be 
written as 0, ll^ 



0.675 



(6) 



where me and mw are the e and W masses. The 
CG interaction cross section is jlij c^^^^p^ — 5.53 x 
10"3^cm2(i;^/GeV)°-3^3. Integrating the flux times 
cross section over an energy bin of say 5.01 PeV - 7.9 PeV 
(which spans 0.1 in logj^Q on either side of the reso- 
nance energy 6.3 PeV), we can calculate the ratio, R = 
^?^/i^f^ + ^tiT'^-, of GR to {v^ + v^,)-GC event rates 
in a water or ice based detector. For an assumed E:^'^ 
neutrino spectrum we get in the optimal case. 



optimal 
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20 



(7) 



In obtaining equation lt7| w e have considered muons with 
contained vertices only |l3 . 

Equation Q, which is a linear function of a;, is dis- 
played in Figures n and m for various values of ^12 (solid 
lines). In these Figures we also show the (dashed) lines 
obtained by allowing small deviations of the angles 6123 
(Fig.nj and (Fig.EJ from their "optimal" values; this 
we do by using the general form of the mixing matrix 
given by Equation Q. 

Figures n and El show that the dependence of the ratio 
R on the mixing angles progressively weakens as x (i.e., 
fraction of origin neutrinos) increases and, in the op- 
timal case, disappears completely for cc = 1 as expected 



FIG. 2: Solid lines: same as in Fig. Q The dashed lines are 
in the same order of Q\2 values as the solid lines — all for 
6*23 = 45°, but 6x3 =9°. 



(see equation 121. Fortunately, for the dominant P7 pro- 
cess (x = 0) the ratio R has significant dependence on 
9i2 thus allowing the possibility of measuring this angle. 
The flip side is that the expected number of GR events 
is smaller in the case of P7 as compared to pp dominance. 

Clearly, the accuracy with which the different mixing 
angles can be probed will be limited by the accuracy 
with which the ratio R can be measured in a real detec- 
tor. While the ratio R always increases with increasing 
012, its dependences on small deviations of the values 
of 023 and 013 from their optimal values are somewhat 
more complicated, though perhaps too small to be mea- 
surable. From the above Figures we see that for a given 
upper limit (say ^ 10%) on the fraction x for a given 
source from independent considerations (e.g., high en- 
ergy 7-ray observations), a measurement of the ratio R 
yields a lower limit on the angle 0i2. Conversely, for a 
given fixed value of 0i2, a measured value of R gives a 
value of the fraction x, thus providing important clues to 
the nature of the neutrino production mechanism in the 
source under consideration. 

At PeV energies neutrinos begin to get absorbed in 
passing through the Earth; this is especially true for Vf. 
at the GR energy where the Ve interaction cross section 
is significantly enhanced. Thus only the downward-going 
to horizontal VeS will be detectable. Fortunately, the 
atmospheric neutrino flux at PeV energies is negligible 
and may not be of concern as a background. The de- 
tectability of Ve GR events has been discussed recently 
for the diffuse neutrino background (due to integrated 
contribution from all sources in the Universe) in Ref. 0| 
for a ICECUBE-type detector. The GR events can be 
identified through the detection of the electromagnetic 
shower produced by the hadronic decay products of the 
W in the detector. This corresponds to ~70% detection 
efficiency |l7j reflecting the branching ratio of hadronic 
decay of the W . The main background, then, are the 
showers due to {v^ + Pe)-CC events within the GR en- 
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ergy window [13 • Interestingly, since, at least in the op- 
timal case, the (v^ + i>e) fraction is equal to the {v^ + v^) 
fraction, and since the CC cross section is essentially the 
same [ig for i/^'s and j/g's, our ratio R in equation {T)) 
multiplied by a factor of '-^ 0.7 gives a good estimate 
of the expected ratio of signal-to-background number of 
events. At the same time, the accuracy with which the 
ratio R can be measured in a detector depends on the ac- 
tual number of GR events detected in a reasonable time 
frame (say 6-7 years), which in turn depends on the flux 
from the source and the size of the detector. 

In this context, we should mention that throughout the 
above discussion we have implicitly assumed individual 
point sources for which the meaning of the pp origin frac- 
tion X is well-defined. For the diffuse background, on the 
other hand, that fraction refers to some kind of weighted 
average over all sources. Our discussions above are then 
valid for diffuse background as well provided the above 
difference in the meaning of x is kept in mind. Clearly, 
from the point of view of signal-to-background ratio, in- 
dividual point sources such as not too distant and yet 
sufficiently powerful GRBs, would be easiest to detect 
because then temporal as well as directional information 
could be used to increase the signal-to-background ratio. 
Also, individual sources are more easily subject to ob- 
servations in other bands such as high energy 7-ray ob- 
servations which may provide independent information 
on the fraction x (assuming hadronic origin of 7-rays 
through 7r° decay). In any case, in view of important 
information regarding neutrino parameters as well as the 
nature of high energy astrophysical sources of neutrinos 
that may be obtained from detection of D,, GR events as 
discussed above, it would certainly be useful in designing 
future high energy neutrino telescopes to optimize them 
for detection of possible i/g's from astrophysical sources 
through the GR channel. 

We thank K. Balaji, R. Mohapatra, and G. Ra- 
jasekaran for pointing out a major error in the motiva- 
tional basis of an earlier version of this paper, which has 
led to a major shift of emphasis in the present version, 
although the main calculations and results remain the 
same. One of us (NG) wishes to thank IIA, Bangalore 
for hospitality. This work is partially supported by a 
NSF US-India cooperative research grant to PB at IIA, 
Bangalore. 



[3] 
[4] 



[5] 



[6] 



[7] 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 



[15] 
[16] 



[17] 
[18] 
[19] 



trophysical neutrino sources for their currently known 
values. 

S.L. Glashow, Phys. Rev. 118, 316 (1960). 
One may also consider the related process iJ^e" W~'y, 
where 7 is a hard photon. Its implications for detection 
of cosmic high energy i^e's have been studied in H. Athar 
and Guey-Lin Lin, Astropart. Phys. 19, 569 (2003). 
See e.g., R. N. Mohapatra and P. B. Pal, Massive Neutri- 
nos in Physics and Astrophysics (Third Edition) (World 
Scientific, Singapore, 2004); C. W. Kim and A. Pevsner, 
Neutrinos in Physics and Astrophysics (Harwood, USA, 
1993). 

For a review and references on astrophysical sources of 
high energy neutrinos see, e.g., I. F. M. Albuquerque, 
J. Lamoureux, and G. F. Smoot, Astrophys. J. Suppl. 
141, 195 (2002). 

P. Bhattacharjee, C. T. Hill, and D. N. Schramm, 
Phys. Rev. Lett. 69, 567 (1992); P. Bhattacharjee and 
G. Sigl, Phys. Rep. 327, 109 (2000). 
J. Ahrens et al. [IceCube Collaboration], As- 
tropart. Phys. 20, 507 (2 004) . 



http 



/ /antares.in2p3.fr71 

/ /www.nestor. org.gr / 1 



http 

http://nemoweb.Ins.infn.it7l 
I. Kravchenko et al. [RICE Collaboration], As- 
tropart. Phys. 20, 195 (2003). 
A. Silvestri et al, astro- ph/041 1007 

\protect\vrule widthOpt\protect\href {http : //www .ps .uci . edu\ 

V. S. Berezinsky, S. V. Bulanov, V. A. Dogiel, 

V. L. Ginzburg, and V. S. Ptuskin, Astrophysics 

of Cosmic Rays (North-Holland, Amsterdam, 1990); 

T. K. Gaisser, F. Halzen, and T. Stanev, Phys. Rept. 

258, 173 (1995). 

F. W. Stecker, Phys. Rev. Lett. 21, 1016 (1968). 
Detection of VeS through the GR channel in ICECUBE 
type detectors has been discussed in Ref. 0|. However, 
they have not considered the possibility of using this to 
measure the angle 9\2- 

L. A. Anchordoqui, H. Go ldberg, F. Halzen, and 
T. J. Weiler, hep-ph/0410003" 

R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, 
Phys. Rev. D 58, 093009 (1998). 

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and 
T. J. Weiler, Phys. Rev. D 68, 093005 (2003). 



[1] S. Eidelman et al. [Particle Data Group], Phys. Lett. 
B592, 1 (2004). 

[2] Unlike the mixing angles, the mass-square differences 
(Am^) appear only in the oscillatory terms in the flavor 
oscillation probabilities which, for cosmologically distant 
sources and for currently known experimental values of 
Am^, average out to zero. The Am^'s cannot therefore 
be probed with high energy neutrinos from distant as- 



